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Abstract: Genetic diversity of species are closely associated with their adaptability and 
potentiality in evolution. Information on genetic diversity and genetic structure of endangered 
species can provide strategies for conservation. In this study, SSR marker was used to evaluate 
genetic diversity and population structure of three endangered Camellia species: C. 


chrysanthoides, C. micrantha and C. parvipetala. The mean number of alleles (3.9), the effective 
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number of alleles (2.328), observed heterozygosity (0.520) and expected heterozygosity (0.501) of 
C. parvipetala were higher than C. chrysanthoides and C. micrantha. In the population level, the 
effective number of alleles ranged from 1.788 to 2.466 and the expected heterozygosity ranged 
from 0.379 to 0.543. The Fs; values between each pair of populations were ranged from 0.143 7 to 
0.453 3, the Nm values ranged from 0.301 5 to 1.488 9. AMOVA suggested that the variation was 
65.72% within populations. Low level of genetic diversity and high levels of genetic 
differentiation were detected for these three Camellia species. The STRUCTURE and PCoA 
analysis showed that 7 populations were clustered into two groups. One group comprised 5 
populations from C. chrysanthoides and C. micrantha, and the other consisted of 2 populations 
from C. parvipetala. In situ and ex situ actions should be urgently implemented as soon as 
possible. 
Key Words: Camellia chrysanthoides, Camellia micrantha, Camellia parvipetala, genetic 
diversity, genetic structure 
1 Biz 
YEN EZR (Camellia chrysanthoides H.T. Chang) (2n230) . /MERHEZS (Camellia 
micrantha S. Y. Liang et Y. C. Zhong) (2n=30) Fl) HSE FEZ (Camellia parvipetala J. Y. Liang 
et Su) (2n=30) EA T RA VU UU FAI = RP TER TEL GK IA AE, 1998; RR 
Mk, 1995). Mee AR T] PEIN ALL, Aes tea T) S] SETH ELA EL 7h 
ME ACAS AY ABP HAL, ETA oo i DX ei HRR. EST ARCRURLR] SP BE TT 
ESB MA, n. VEIN AA ne LE MEL, ABODE Ee SS eH AE, ix 
3 “MAP EE ETEK RE IR, FRE MO ERE RTETRRUIVE EA CC CP 
Ta EL 52 RO A se) ABH eh ESE TS, 2017). 
3 f S PEPE Re a Rb ET ET, EE IW PS eA ee (Frankham et al, 
2002). ED Pr PH Se ELD TE A PEAT TA PT FE ESS BO BE OLA, 3S 
AY DLA DP SRS 55; EEN LIT itil] SE HE DERE Te FEI (Segarra-Moragues et al, 2005; Su et al, 
2017) IRI T b iG A Srt OE PEE A ZT OAR YO IE SF 2 125 (Ryall, 
1998). $4. P/£ CMicrosatellites) pric th Hl] fa] ERFA] Csimple sequence repeats, SSRs) , 
Hom ed. JOE. EMERE ERRER EI, EMT 
FEREZ- UR ICE RA Li et al, 2013; Gyorgy et al, 2014; Meng et al, 2015). $2 BA id) 
i£ MB T So REUS S REPERI FE (Yang et al, 2016; Ma et al, 2015; Zhao et 
al, 2017; He et al, 2017). 
RS TED He WA BUG an, "E e ABE ERI A I UR 
HEHE HE BUR lS (Cires et al, 2011). FPR EAE. “Tes TEAR AT SIE Se TEAR TG E FEE 
FUEL OT FAS DLR» AE, ASE FEAF (ICE bi th = EATER 7 T RERBA S 
PEPE MURAI, ETE SARE TAG EREEREER, ART TAR de D PRP SES 
2 MRI 
2.1 HR 
TEE UE diz 46 AR M E Sz 48 AR TE 2) KERA D BC GRE IK A AE SE HH, 1998; Ming & 
med 2007), Ming & Bartholomew (2007) 7) JE 446 25 AIF BUTE TEAS (Camellia 
indochinensis Merrill) , [HZé/] ES TEAS UAT E AER TETAS LFE EKA Kk, AS 
FETE ATK AR IK AE HE (1998) TEL St Se AE AS VE IPSE SEE AT RE RITENA 
ARRE 7 AARE, 3 DDE TEASE 2 T SIE Se EAS EA 2-4 RT AS RE Jo 
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tats. FREE BHI fU BH2 WESINTEN ELESJSCK HU. 3cm, 5 CR Bil GKRIAAME 
HEU, 1998)J8XR f EMT AETEZR TEM BEA 4-5.5cm AFA. [Hil ETE ACE 
R, NURDAECARN DUCH ACBUTCGIKCEXA, 1979), 46 HUE 4-5.5cm JEU UE EET Hi 
ROHR AREAS, 2013); BH2 PEALE PA T 4628 ZA V d ES PE B SEE RAE A S| b 
FE BHI Rn HEE SERE EATER IR SE EI A LL, RBA A SE Ix PS PEE 
AVE ETERS © EPR YI BX 2-3 Jr Br Hee ARH TBA AE ERE E BE SS Ra. HAFA 
DNA WIPER. PRIUS. SEED AME SEL 1, BPR ALA 1. FEUER ANE) 
FEW FLT ATH (BK) . 
#1 —BRREGUCRTAORECKSE ag A 


Table 1 Sampling information for 7 populations 


AR FPFE 
VURREuERRAR D) P Longitude ET Ky) 5 — 
Species (Voucher Popul KEHA Location (Œ) ^ Altitude Popu ee 
number) ation Latitude (m) ation 
(N) size 
EHEER BH1 EIEI E BN 106?44"/ 305 «100 34 
C. chrysanthoides Anzhen, Bingiao, Longzhou 22?15' 
(2017042201) BH2 EME ERE S BEIT 106°47'/ 251 <50 32 
Zicong, Shangjiang, 22°10' 
Longzhou 
ieste A NB1 BETRA BE ASI 106°51'/ 337 8 8 
C. micrantha Fengle, Xiashi, Pingxiang 22°09' 
(2014011601) NB2 A ETE 2 AERA 106°52'/ 370 12 12 
Fengle, Xiashi, Pingxiang 22°09' 
NB3 Sere A AY 106°51'/ 355 «100 35 
Xiashi, Xiashi, Pingxiang 22°07' 
ESE ALAS PT1 THY BURR 2 IR GY 107°08/ 208 <50 34 
C. parvipetala Paitai, Zhilang, Ningming 21°50' 
(2014011801, PT2 TAH ECTS EC A] 107213"/ 492 «50 29 
2014011802) Name, Tongmian, Ningming 21°49' 
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Fig.1 Map of the sampling location of 7 populations 

2.2 DNA FERMI SSR 2) 7 

KHAR AY CTAB (Doyle, 1987) fe Ut Fr A DNA. M Lu et al (2014) fl Liufu et al (2014) 
ZR AEAEAS PUR BE BAER T ASE 59 o6 pd Dg S| arp aE 1106] D RST, 2 ASHE 
ral) CELER S| OA AF. PCR PMA RAPE 75 Liufu et al (2014) E AMT. 
2.3 BE HT 

Git BELA ed BHAA Genepop version 4.1(Rousset, 2008) fall) Her itt -YAI ^T^ ft 

(Hardy-Weinberg equilibrium, HWE) , Jf] Prf$ P [E XtfT sequential Bonferroni f IE(Rice, 

1989). GenALEx 6.5 (Peakall & Smouse, 2012) X fF HF Stitt -F JS ie A CN) . ARSE 
PFE AA (Ne) . WUMARAEE CHo) . WERE He), ERR (OO JUGE 
HD (PPB) SEA REVETERL. A 1 28871224) 9T AMOVA (Analysis of molecular variance) 
(Excoffier et al, 19928] H] Arlequin 3.0(Excoffier et al, 2005) f IFTE SE, HAH APE ER VER A] A 
TREE ASI Fer. SEA (Nm) fh SER AA M = (1 - Fsr) /4 Fsro 

EG FEE AE LA) aT CPCoAO IH] GenALEx 6.5 SX fFi- T£. STRUCTURE 2.3 
(Pritchard et al, 2000) %t 7 “S#EREUEAT 2191, TR DEEP MOT BU 2E REO TPA. RI 
WEAN: K=1- 6, ^A K 8^) 9i fT 20 MW, Burn-in 10° Xx, Markov Chain Monte Carlo(«MCMC) 
5X10 IR. ISAT SERA AA ERG HE STRUCTURE HARVESTER (Earl & Vonholdt, 2012)) 
Br vU ERU FE ZAS K. FAY GenALEx 6.5 FF) Mantel test W 7 APPRENDE EEE A 
EBSA Cisolation by distance, IBD) FEIKE. 
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3.1 fiiius e ETE 

11 ML ES RAE AR 184 4 HH UI SS AS br A [A] 92 4S, E35) REA hr 8.364 NR 
2). T4 TREE 11 in SE 77 1X Hardy-Weinberg ^F ES 3$, Eds JE 9 Uefa a (P 
< 0.05) , 2 sequential Bonferroni EJ, RAMA TERS ÆRE NB3 (iAP. PUR dS 
A FA Ja BEd iF 


X 2.11 ML eoe S PPE a S. 
Table 2 Population genetic parameters of 11 microsatellite markers 


WW ap ec ae | —— 

ee il WMA WERA RHR ERM I HE LA E 
Loci N, N, N; E Ho  EHe I Fis AR Fit RRM Fsr 
FLA2 9 39 2283 0.592 0.535 0974 -0.107 0212 0288 


FLA7 9 4.7 2.635 0.605 0.602 1.118 -0.005 0.206 0.209 
FLA20 12 4.3 2.741 0.591 0.567 1.060 -0.042 0.296 0.325 


FLA23 9 5.0 2.492 0.626 0.567 1.092 -0.103 0.231 0.303 
FLA27 6 3.4 1.980 0.436 0.435 0.757 -0.002 0.372 0.373 
FLA30 5 2.1 1.653 0.286 0.279 0.457 -0.024 0.302 0.318 
FLA32 10 4.3 2.157 0.521 0.495 0.904 -0.053 0.144 0.187 
TER3 5 2.7 1.541 0.306 0.315 0.546 0.028 0.298 0.277 
TER8 7 2.7 1.777] 0.320 0.366 0.624 0.126 0.481 0.406 
TERIS 12 5.3 2.428 | 0.643 0.567 1.108 -0.135 0.100 0.208 
TER21 8 3.3 2.083 0.430 0.448 0.786 0.040 0.407 0.383 


“FY — 8.364 3.8 2.161 0.487 0.471 0.857 -0.025 0.277 0.298 


Mean 


TE: Nr ASME AR Ny. FIA EA, Ne AAEE Ho. WMA EE; He. HUA AE; 
L ARM: Fis. VAR ZR Fit. APIA A BM: For REAL IRL 


Note: Ny means total number of alleles; V4 means number of different alleles; Ng means number of effective alleles; 


Ho means observed heterozygosity; He means expected heterozygosity; / means Shannon's information index; Fis 
means inbreeding coefficient; Fit means genetic differences among entire population; F;; means genetic 
differences among subpopulation. 

TORRES & FPL TA BUT EAS Le 3. EMUIOE E, “ails HEAP EJ Sr e SUR 
(Na) « ARFER ON). WYRE Ho), BAR AE He) IIE 3.9. 2.328. 
0.520. 0.501, TT TETTE SETEZR RU] Tee TE AS | EREKE, 7 ARIES EE DIC CN) 
H 2.7 (NB1) ~5.1 (NB3) ZLIR], PHIEN 3.8, ANEMIA (Ne) TE 1.788 (NBI) ~ 
2.466 (PTI) Zia], “PHYA 2.161. Willa Gee (Ho) Æ 0.409 (NBI) ~0.543 (NB3) 
ZH, PHEN 0.487. WAR AE (He) 1E 0.379—0.543 ZE, “PHBA 0.471, REE 
A NB3, AUÉKTERREENBI (43) - 

de3 SPEEA 7 NEE ETE K Z 


Table 3 Genetic diversity of 7 populations 


Vz ^ ds vm 

«0013. qu X^ un 200 50 

Vh Species RES Stl ey pe ma PUERO MBA 
d Population 3&4 ES ae — GV OF P 


Ht N. Ho EH 
HN, ANa Eio. EHe PPB(%) 


EU AE AES BHI 3.1 1.914 0441 0.411 -0.066 81.82% 
C. chrysanthoides BH2 4.2 2.064 0.452 0451 0.013 90.91% 
EY) Mean 3.7 1.989 0.446 0431 -0.024 
INERTES NB1 2.7 1.788 | 0409 0.379 -0.114 81.82% 
C. micrantha NB2 3.5 2.349 0.523 0.508 -0.041 100.0096 
NB3 5.1 2.356 | 0.543 0.543 0.027 100.0096 
IJ Mean — 3.8 2.165 0.492 0.476  -0.038 
1) HME BE AE AS PT1 4.3 2.466 0.508 0.540 0.063 100.00% 
C. parvipetala PT2 3.5 2.190 0.533 0.463 -0.139 81.82% 
I Mean 3.9 2.328 | 0.520 0.501 -0.028 
^35 Mean 3.8 2.161 | 0.487. 0.4471 -0.037 90.91% 
Es Ny. SRR, Ne tC EDLAEDUÉG Ho. WUWZRERE; He. WARA: F. DEUS PPB. 


AS ABA. 


Note: N4 means number of alleles; Ng means number of effective alleles; Ho means observed heterozygosity; He 


means expected heterozygosity; F means fixation index; PPB means the percentage of polymorphic loci. 
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65.72% (44) . TERRI, CNWD TERRA UA RETE 


AE IR] 


SEPIUS, EIERE EE i 24.4190, BREA RSH 75.59%; “Fee de ZR PIE IR] AE 


AJ 


18.2096 , FEE A AERE E 


t1 


81.8096; “HE TEAS PE 


IARE i 35.39%, HAE A AEH 64.6196. 


V Pa Ph EEE) IE LAR Psp MUZE AT Nm WZ 5 AAS, Fsr Æ 0.143 7~0.453 3 ZA), HE 


HA, 1 ZEURPREE RI Fore) F 0.15; 6 H7 
T 025, DLAA RRR. Nm 1E 0. 


301 5~1.4889 Z [HR], (WA 3 ZEUPEREETR] 


HEEJ Fr (HE 0.15—0.25 ŻE), 14 ZEUREREIR] Fer EC 


ZET 


ES RMAF 1, 


PREE AMER. TEM se TE ARAM AE Se EAR RI A OD RL F sr 73 0.173 1, HEERA 


PRESALE ARATE LAB Fsr 73 0.258 3, 
TRU TERE EER ALB), De TEARS 


HA pA 


R4 DPR CAMOVA) 4558 
Table 4 Analysis of molecular variance (AMOVA) for three Camellia species 


FE SEALE ARAN IE TEAS 2 15 83. F sr 73 0.206 8, 
BAER A MLBEK 


RFRA EE EI 


4y 2H. Grouping SERA . Percentage of {LAM 
Source of variation i» 
variation (96) Fsr 

BUE AMA tla] Among species 9.66 0.34** 

All individuals Th AY PREET] Among populations 24.62 
within species 

PPIE Within populations 65.72 

EHEER TiEIR] Among populations 2441 0.24** 
C. chrysanthoides PPIE Within populations 75.59 

\ FEBAERE FPEF] Among populations 18.20 0.18** 
C. micrantha #t#E A Within populations 81.80 

SHEE TEAS FPEF] Among populations 35.39 0.35** 
C. parvipetala PPIE Within populations 64.61 


ik: "XR P «0.01. 


Note: ** Means P « 0.01. 


d 5 TUBERISUPEAMESRÉ For OG fi PO AIDE Nm WARE) 


Table 5 Pairwise comparisons of F7 (below the diagonal) and Nm (above the diagonal) for 7 


Camellia populations 


BHI BH2 NBI NB2 NB3 PTI PT2 
BHI — 0.774 1 0.399 4 0.611 1 0.533 7 0.402 7 0.301 5 
BH2 0.244 1** — 0.887 8 1.260 0 0.761 2 0.562 6 0.333 9 
NBI 0.385 0** 0.219 7** — 1.489 9 0.728 4 0.580 2 0.381 5 
NB2 0.290 3** 0.165 6** . 0.143 7** — 1.301 2 0.791 4 0.439 3 
NB3 0.319 O** 0.247 2** | 0.255 5** 0.161 2** — 0.586 8 0.375 1 
PTI 0.383 O** 0.307 7** . 0.3012** | 0.240 1** 0.298 7** — 0.461 4 


PT2 0.453 3** 0.428 1** | 0.3959** . 0.362 7** ^ 0.3999** 0.351 4** — 
YE: "RAN P«001. 
Note: ** Means P « 0.01. 

STRUCTUREZIPCOoA 4] jr = fi ae HE ASG JB S REAR STRUCTURES 2 
RER, TAMREISAA SAT ae EEE) K=2 (E20 , BEI PT See RO PAL) 
TERETE ASST AREER EA DB DA — 28; ZEAE IST PIRE B Aad PRA A — 28 R13) o 
PCoAS Pr SR CIAO , TARARE AT PRA APA ZE, BUY HESS 2 Sh REA) Se FES 
HJ34  FEUEEA A ZA, ZI TEAR -S RE) 28, (EL I TER EH SL 
P, RRREPTISE RET FN FES TERRA S RE. Coord. 1 FlCoord.24) 3&3 18.76% FI 12.59% HY 
AER. Mantel test 22284 CRIS) , TA REET SREB PS A ee FB ES E 5 ASE AHR PEE 
iE (R7=0.0847, P=0.23) . 


Delta = mean(|L"(x)]) / sd(L(x)) 
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Fig.2 Relations between the number of determined group K and estimated value AK 
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Fig.3 Two clusters were detected from STRUCTURE analysis. 


= p ^3 235i Principal Coordinates (PCoA) 
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Kl4 -BRARACARIATA WI (PCoA) fi 
Fig.4 Principal coordinates analysis (PCoA) of SSR phenotype of three Camellia species 
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Fig.5 The correlation between genetic distance and geographic distance for 7 populations 
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4.1 RETETE 

ip E EVE AE E E PE IE TS J AY E AA (Frankham et al, 2002). — R2 TG, 
KATE. BAIE ARIO BUE RAS RURKI E & FE E(Spielman et al, 
2004; Gao et al, 2017). -5 i] f AH EG, SMEER FEMEA N MERSE He 
BUR GEHIE NAT 3.7, He =0.431; MERER, Ní4 23.8. He 2 0476; NIERE, 
N423.9, He=0.501) » {RFR WIRE ER CC. flavida, A CSEAISEDNID =4.4, He 
= 0.555) (F X4, 2015), KEEFE (C. taliensis, Ag CEMER BE) = 6.776, Hs JEANS 
FETE) = 0.597) (Zhao et al, 2014), EFI CC. japonica, A= 16.5, He = 0.84) (Ueno et 
al, 2000)#ZE (C. sinensis, A=4.3, He =0.64) (Yao et al, 2012). FAL, =A EASY I 
BAR BOR ACE B FE FETE. ARRE PEPER RAS FH, PG: AE APIS, HE 
4) TRUST] SX (Nybom, 2004). — fia TEA 2) THAR Fs RE, FF AAA A TER 
HEEMI HE AAA, SBE TE RE) BME, SERS ha, MARRE 
SE FRET ABD F 100. f HABE. PHREN Fr Tad ti RT SB SS EE 
WES FETEZCOTHIE o 
4.2 RAE 

= Fh ae eA PRE FZ CER CE 5S), EE, EOL ARR Pop MEE 
NB1 4I NB2 ŽEF 0.15, FERSIERE BOR FE. EZR TEREX C0.15«6 H 
PREE For 18«0.250. EXTR KR EE2M (14 PREE] Fs; 1850.25) (Wright, 1968). FFA 
PRE AI EVE ETERS (0.244 D) FFE FEAR (0.143 7, 0.255 5 800.161 20. ERRERA 
4H. 7) URE TEAR AT A PP RE TE) ARCA B3 HA, (0.351 4) , AFAR ERA PY RE AE 
AYE, FARE AD IR TEASE PT2 AL BHR BESTE AN AREER. PARK AE 
FE ASAE FRE Fd Ja CRUEL AE P Md BC AK, 2015). PRR AEA 
HURAE, Ford ps in e RR PEER A EL SR — (H7) 83, 2000). 7 AEREA 
Alii Nm f^] (425) , BUB 3 APPR Nn AF 1, AG Wright(1931) i6, FUE HPR 
[RJ Nm»1 BE, ERRA Be IER FERAE RITE FI SEB LE PS Yh PE RI XE P OP BS A 
^E. FERN OR BE Pe HR AE BE SK. FERRO, AG 
{EAR (Wei et al, 2017; XENARJI, 2017). WHACKMEAS, 2010) MAHA (Liu et al, 
= 2012 ANH FEE RETR ED AIR, SBC BRA BA UP SEAL, HERRN. — 
oa PRR TEAS EE Sr ha. BR TE AEE PT E SSPE TB] FE BER E 
© ZEA, ETE PRR) ADIL, See MBE K . PRE Ta ARAE 

FRAGA. PRR PUG EER ed, YE SN Se EAS K ET a fe AME. 

STRUCTURE $I PCoA DIRKA, URE AEE a fe op A 2, Bos ae te 
ARAN FE TEAR AEB MMA AN 2A, "IEEE KBAR A, Se — se PEAS 
BL} Ati LX Sa $86, A ERE RATERS I ETE AR ZAI BU), IER TEARS FCR PE EAE 
ASA MOBEK , Be BH EI EAE AR B 2 (SEDO HARE AY AD TE LER E] 3 SRR HEARS RT Be e 8] — FP 
HY 
4.3 tPA SL 

YDB AR E FEVEZKOT E RET Be 0 Pa NE 7J 95 YA (Hamrick & Godt, 1996). AWE FE 
SERA, = RIA FOR RRK EER & EE ER OP A RP BPP ed Er 
3, PEPR E Sa eB COS BAS Ta Ee, ARE PT2 opi Hh SOT AC REENA 
Bt; FRE BHI A 7h CAROT RARER; FREE BH2 DREA, ORT 
Sk. ARF BH2 EZ, SARKARI. PAK, MAA PRE 
4f Sic a fe US PR HE RP BE HRP 6, SKE DISP TREE RE RP AB ER 
d EA PRXENCRIULVEUS Pe RP A 
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